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Entangled multi-qubit states may be generated through a dispersive collective QND measurement 
of superconducting qubits coupled to a microwave transmission line resonator. Using the quantum 
trajectory approach, we analyze the stochastic measurement traces that would be observed in ex- 
periments. We illustrate the synthesis of three-qubit W- and GHZ-states, and we analyze how the 
fidelity and the entanglement evolve in time during the measurement. We discuss the influence of 
decoherence and relaxation, as well as of imperfect control over experimental parameters. We show 
that the desired states can be generated on timescales much faster than the qubit decoherence rates. 



INTRODUCTION 

The realization of quantum- optical concepts in con- 
densed matter systems has led to remarkable progress 
during the past few years. One of the prime examples 
is the study of quantum electrodynamics (QED) in su- 
perconducting circuits. Earlier suggestions to implement 
the Jaynes-Cummings model in the solid state [U EJ |3j 
were followed by a proposal [4] to employ on-chip mi- 
crowave resonators and couple them to artifical atoms in 
the form of superconducting qubits. This seminal idea 
was soon thereafter realized experimentally [5|, creating 
a solid-state analogue of conventional optical cavity QED 
[6] . The tight confinement of the field mode and the large 
electric dipole moment of the "atom" yield extraordinary 
coupling strengths. As a result, these highly tunable 
systems have been employed to demonstrate experimen- 
tally a variety of achievements, including: The Jaynes- 
Cummings model in the strong-coupling regime [U |7l |8] , 
Rabi and Ramsey oscillations and dispersive qubit read- 
out O [To], generation of single photons fTP and Fock 
states [m [13], cavity-mediated coupling of two qubits 
[HI [15], setups with three qubits Berry's phase [17] , 
and the measurement of the photon number distribution 

The strong coupling makes dispersive quantum non- 
demolition (QND) readout possible, both for qubit states 
and for detecting single photons [19]. QND measure- 
ments are ideal projective measurements that reproduce 
their outcome when repeated [20l [21]. Any QND mea- 
surement may be applied to (probabilistically) generate 
states. In particular, having several qubits inside a com- 
mon cavity (as realized in recent circuit QED experi- 
ments [HITS], for a schematic setup see Fig. ([l]))), one 
may produce entangled multi-qubit states, even without 
employing directly any qubit-qubit coupling. In the con- 
text of circuit QED, this option has been investigated 
previously in a series of remarkable studies [22l [23l [24] . 
However, these consider primarily two qubits, with a re- 
cent work [23 discussing the extension to more qubits 
in general terms. The present paper aims to go beyond 
these studies in several aspects. First, we present nec- 



essary conditions for being able to generate arbitrary 
multi-qubit states out of a given subspace of the total 
multi-qubit Hilbert space, using only single qubit oper- 
ations and subsequent collective measurement. Second, 
we carry out detailed quantum jump trajectory simu- 
lations also for the case of three qubits, where W and 
GHZ states may be produced. We show how entangle- 
ment is generated in the course of the measurement pro- 
cess, paying attention to the effects of relaxation and 
decoherence. Moreover, we analyze how imprecise fine- 
tuning of experimental parameters would lead to a loss 
of entanglement after its initial, transient generation. Fi- 
nally, we comment on possible experimental realizations. 
Such a measurement-based scheme complements other 
approaches for entanglement-generation in circuit QED 
[l|23|2i|27l[28l[2l[3a[nil32] , based on unitary dy- 
namics, and may prove advantageous for some purposes, 
since generation and measurement are combined into one 
step. It might also be used to generate entanglement be- 
tween qubits in spatially separated cavities, without any 
direct interaction. 

MODEL 

We investigate a QND scheme utilizing the coupling 
of superconducting qubits to a bosonic field mode of a 
microwave resonator as examined in [4l[5].The presence 
of excitations in the qubits inside the cavity gives rise 




Figure 1: (Color online) Schematic setup: Three supercon- 
ducting qubits (red) are coupled to a mode of a coplanar mi- 
crowave resonator (blue) . The measurement of the phase shift 
of a transmitted microwave beam can be utilized to rapidly 
synthesize e.g. maximally entangled multi-qubit states like 
GHZ and W-states. 
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to a frequency shift of the cavity mode, which can be 
observed dispersively via the phase shift of a transmitted 
beam. In turn, the measurement backaction leads to a 
projection of the qubits on a state that depends (a) on 
the chosen set of couphngs and (b) the initial (product) 
state the qubits are prepared in. 

We consider a system of a driven cavity mode coupled 
to N qubits 
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a model commonly known as the Tavis-Cummings model 
which has been recently realised experimentally for = 
2 [HI [151 and TV = 3 [16 . The first term of this Hamil- 
tonian describes the cavity mode with a frequency cj, 
the second all qubit energies, the third term realizes the 
Jaynes-Cummings coupling for each qubit to the cavity 
with bare coupling constants while the last term de- 
scribes the driving of the cavity with the readout mi- 
crowave tone which will yield \e\^ photons in the res- 
onator on average (/^cavity is the intensity decay rate for 
the cavity). 

In the limit where all the qubits are strongly detuned 
from the cavity, it is well-known [4 that the qubits impart 
a state-dependent phase shift on the cavity mode and the 
effective Hamiltonian can be written as: 
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energies could always be chosen very different, such that 
this unwanted interaction does not play a role, since the 
qubits are non-resonant, (ii) As we will argue later in 
more detail, the measurement rate V is usually much 
larger than thus making the effects of the inter- 

action negligibly small even when the qubits are in reso- 
nance with each other, (iii) We note that all the states 
we consider as examples are eigenstates of the flip-flop 
interaction, Eq. ([3|. Therefore, even if the qubits are 
chosen to be in resonance (as is ultimately assumed in 
our simulations), the interaction will not have any im- 
portant effect on the dynamics besides trivial phases be- 
tween subspaces that are rendered mutually incoherent 
by the measurement anyway. Thus, we will neglect the 
flip-flop interaction. 



NECESSARY CONDITIONS FOR THE 
GENERATION OF GIVEN TARGET STATES 

The protocol we are envisaging is to first produce a 
product state, using local operations on the individual 
qubits, and then to project on an entangled state by mea- 
surement. This scheme will be successful with a certain 
probability. Our aim in the present section is to briefly 
discuss the necessary conditions that must be met to be 
able to generate a given class of entangled states. It goes 
without saying that once an entangled state has been 
reached, one may then apply further local operations to 
reach a corresponding subspace of the full multi-qubit 
Hilbert space. 

Let us first fix notation. The coupling strengths Qi de- 
termine the phase shifts induced by the individual qubits. 
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and for convenience we will collect them into the coupling 
vector 
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The Hamiltonian ([T]) also induces an effective flip- 
flop interaction [4p ] of strength J^^^ = gaSfdi^a + 
A/3)/(2AckA/3) between each pair of qubits in the 

same cavity (for couplings ga{i3) and detunings from the 
cavity Aq;(/3), in the dispersive limit \g\ <C |A|): 



where the overall strength g just determines the measure- 
ment time-scale but does not affect the reachable states. 

Using qubit excitation operators hi = ^^^^"""^ , we de- 
fine the measurement operator TV as 
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When simulating the master equation (11) to be derived 
from the Hamiltonian ([2|, we neglect this interaction for 
several reasons, (i) In a concrete experiment the qubit 
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Note that even in current experiments, the gi are tunable 
in magnitude and sign simply by choosing the detuning 
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A appropriately. The frequency shift imposed on the 
microwave cavity wih then be 
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where we have defined Ui := (hi) G [0, 1] as the excitation 
number of the i-th qubit. 

The desired entangled state |^) = has 
complex amplitudes 
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in the energy eigenbasis of the qubits (a product ba- 
sis that diagonalizes fii). Thus, we have rii^^) = 

In the following we derive necessary conditions for be- 
ing able to produce arbitrarily chosen states out of some 
M-dimensional Hilbert space that is spanned by a sub- 
set of M basis states In order to generate a certain 
target state given by arbitrary c?, we need to adjust the 
couplings such that all base kets with non-vanishing aj 
yield the same phase shift. Assume the amplitude vector 
of the target state has M G {!,..., 2^} non-zero entries 
aj^j G {1, .., 2^} where the corresponding indices can be 
written as a family Fq, with dim(Fc^) = M. Then the goal 
is to use the measurement to project the system onto the 
subspace given by span ({l^/^'j) |j G Fq, }). 

In the simplest case this is directly possible by choosing 
(i) an appropriate initial product state of the qubits (to 
fix the amplitudes) and (ii) a suitable coupling vector (to 
project onto the correct subspace). 

Choosing an arbitrary initial product state allows for 
the choice of 2N complex amplitudes. Due to normal- 
ization of the TV single qubit states and a an arbitrary 
global phase for each of those states, we essentially have 
2N real parameters to choose. 

The amplitude vector of the target state will - up to 
a constant common factor due to the renormalization af- 
ter projection - be determined by the amplitudes of this 
initial state. This suggests that, in general (i.e. for arbi- 
trary target states), we can only aim at reaching states 
that satisfy 
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Again, we had to subtract 2 to account for the irrelevant 
global phase and normalization. 

Note that for the maximal value of M = 2^, the last 
inequality does not hold for any > 1 and we recover 
the fact that arbitrary states are in general not prod- 
uct states. Note that we have just found a necessary 
condition for constructing arbitrary states out of an M- 
dimensional subspace. When choosing particular states, 
e.g. trivially separable states, one may still be able to 
construct those even if they formally violate Eq. Q). 



We now turn to the question when it is possible to 
choose the couplings such that the measurement cannot 
distinguish the components of the target state from each 
other. This requirement of equal phase-shifts formally 
corresponds to a set of M — 1 equations 



(10) 



where i,j denote successive indices out of Fq,. 

As tunable parameters to our disposal we effectively 
have A' — 1 couplings (discounting the overall strength g) 
so this set of equations will in general be solvable as long 
as M < A" is fulfilled. 

As we will demonstrate below in several examples, 
some of the most interesting entangled states, such as 
W and GHZ states for three qubits can be synthesized 
by this scheme. Indeed, they have M = A" for the W- and 
M = 2 for the GHZ-states and thus satisfy the necessary 
conditions discussed in this section. 



STOCHASTIC MASTER EQUATION 

In this section we turn to the quantum trajectory ap- 
proach known from quantum optics [21, |33[ |34l 135] [36l 
EZllMlEiliOlliTll^liaSl]. The stochastic master 
equation to be presented below allows us to model the 
backaction of the phase shift measurement and to pro- 
duce individual realizations of the measured phase shift 
signal. 

In the case of a cavity mode that decays much faster 
both than the qubit decoherence rates cavity ^71,70) 
and the couplings to the cavity mode tvcavity ^ ^iVi G 
{0,1,.., A"}, it is possible to adiabatically eliminate the 
cavity mode from the system and find for the stochastic 
master equation (in the interaction picture) for the qubits 
alone 
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Here T 



is the measurement rate, 71 and 

are the qubit relaxation and dephasing rates. Pi is the 
projector onto the excited state of qubit i, and \e\'^ is 
the average photon number circulating inside the cavity 
mode. See [42l |45] for a detailed derivation, and [46 for 
our recent analysis of photon detection in circuit QED 
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using the same approach. The stochastic master equation 
is conditioned on the measured signal 
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where ^ represents the fundamental, unavoidable vacuum 
noise (with {^{t)^{f)) = 5{t - t')). Physically, X{t) is 
the appropriate (suitably normalized) quadrature com- 
ponent of the electric field transmitted through the cav- 
ity, which is proportional to the phase shift that indicates 
the multi-qubit state. Experimentally, this signal would 
be measured in a homodyne detection scheme. Note that, 
for a two-sided cavity, information is contained both in 
the transmitted and the refiected signal, and we have 
assumed that both parts of the signals are superimposed 
symmetrically to extract the maximum possible informa- 
tion content [46j. 



EXAMPLES OF MEASUREMENT-GENERATED 
ENTANGLED STATES 

In this section we discuss the most relevant examples 
for the case of two and three qubits in the cavity. More 
precisely we will show that it is possible to generate Bell 
states, W states, and Greenberger-Horne-Zeilinger states 
(GHZ-states). 

We will be able to observe that the measurement in- 
deed first drives the system to one of its attractor solu- 
tions (among them the desired state) which are then sta- 
bilized by the measurement. The attractor nature of the 
subspaces selected by the coupling vector ^ can be im- 
mediately understood from the structure of the stochas- 



tic master equation (11), by realizing that the stochastic 



term and the measurement induced dephasing term both 
vanish if the density matrix is in the desired state. Only 
relaxation and dephasing can take the system out of this 
final state, and we will discuss their infiuence later. 



Quantitative characterization 

In order to characterize the time-evolution during the 
measurement process, we have plotted several quantities. 
We plot the phase shift signal X(t) and the excitation 
number in each qubit as functions of time. To verify 
that we have indeed obtained the desired state, we will 
compute the state fidelity F between the density matrix 
from the simulation, psimi and the ideal state density 
matrix a, according to F = i'c\■^fp^j^a^^J~p^i^\. Finally, 
the two-qubit entanglement between two qubits A and B 
will be measured by the log-negativity. Given the density 
matrix p of the two qubits (after tracing out other qubits, 
if needed), this is defined as En{p) = log2 Here 
p-^^ is the partial transpose with respect to qubit A, and 
\\p\\ = tr[^ p^ p] is the trace norm. 



Before discussing the individual examples, we briefiy 
point out the general features. Looking at the results 
[e.g. in Fig. ([2|], we find that in those cases where we 
end up in the right state the fidelity as well as the log- 
negativity are 1. Furthermore, the state is stabilized by 
the measurement, meaning that, due to the absence of 
any non-vanishing terms in the master equation's right 
hand side, it is frozen. We observe that the state is gener- 
ated on a timescale given by the measurement rate F"^. 

While discussing the examples we will also analyze 
plots that show the probability density of various quan- 
tities evolving over time. This point merits a brief dis- 
cussion. The time-evolution of the distribution for any 
simple quantum-mechanical observable can be immedi- 
ately obtained from the time-evolution of the average 
density-matrix, i.e. from the standard, non-stochastic 
master equation. In that case, simulating a large number 
of stochastic trajectories and then averaging over the re- 
sults would be unnecessarily cumbersome. However, that 
argument becomes void as soon as one considers signals 
that depend on the entire pre-history of the trajectory. 
An important example is the time-averaged cumulative 
phase-shift signal. 
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This quantity has the advantage of tending towards a 
well-defined limit in the course of a QND measurement, 
with the fiuctuations around that limiting value decreas- 
ing like 1 / • It is not possible to obtain the distribution 
of X from the average density-matrix p, and quantum 
jump trajectory simulations are needed. 

Another example is represented by quantities that de- 
pend non-linearly on the density matrix. In those cases, 
the average density matrix is irrelevant since, obviously 
(/(p)) 7^ f{{p)) for a nonlinear function /. An impor- 
tant case is the entanglement measure En- In fact, the 
average density matrix is never entangled {En{{p)) = 0) 
for our examples. Thus, it is indeed necessary to obtain 
£^ AT for a large number of trajectories in order to discuss 
its statistical behaviour and plot the probability density. 



Bell states for two qubits - no decoherence 

In the case of two qubits and vanishing decoherence 
rates 71,70 = 0, the generation of Bell states is straight- 
forward [23 . We imagine starting the experiment with all 
qubits in the ground state |00) and applying a Hadamard 
gate {7r/2 - ax rotation) at some time to, which leaves the 
system in the product state |^o) = FI^ + l^)) — 

\ (|00) + |01) + |10) + |11)). We want to generate the 
Bell state 
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(14) 



which is the two-qubit version of a W- state. Clearly 
the amplitude vector for this state is simply V^c? = 
(0, 1, 1,0)-^, and the resulting Eq. (10) for the couplings 
is given by gi = ^2, thus ^ = (1, 1)-^. The desired state 
will be generated with a success rate rj given by 



T] : 



(15) 



meaning that the experiment will in 50% of all runs end 
up in the correct state (as confirmed by observation of 
the correct phase shift). 

Likewise, for the Bell state vector |$+) = 
^ (1 00) + 1 11)), we find for the amplitude vector >/2 = 

(1,0,0,1)-^, and for the characteristic equation for the 
couplings gi = —^2 ,which is fulfilled by the choice of cou- 
pling vector G = (1,-1)-^. Note that in principle \^~^) 
could also be generated by first producing |^+) and then 
applying local unitary operations, and the same is true 
for the two other Bell states, |^~) and |^~). 

Individual traces and probability density time evolu- 
tions for various quantities are shown in Figs. |2]and|3) 
respectively, for the two types of Bell states discussed 
here. 
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Three qubits - no decoherence 

Generation of W-states 

Similarly, for three qubits, the generation of W-states 
is straightforward as well. We imagine starting the ex- 
periment with all qubits in the ground state |0), applying 
a Hadamard gate (7r/2 - ax rotation) at some time to, 
leaving the system in the state 
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|010) + |011) 
|110) + |111)] 



We aim to generate a W-state which for three qubits is 
given by 

(1001) + 1010) + 1100)). 

We find the corresponding amplitude vector \^~ct = 
(0, 1, 1, 0, 1, 0, 0, 0)-^, and the resulting equations for the 
couplings, gi = g2 = 93^ solved by equal couplings to all 



Figure 2: (Color online) Generation of the Bell state |^^) = 
(|10) + |01))/y2 (the two-qubit W-state): (a) Quantum tra- 
jectories illustrating the different phase shift signal traces 
X(t). Three traces have been selected, corresponding to the 
possible outcomes of the measurement given the same input 
state 1^0) • At time Tto = 3, Hadamard gates are applied 
to both qubits, starting from the ground state. As in ev- 
ery real measurement of field quadratures, the signal X{t) is 
smoothed by doing a windowed average over a suitable time- 
span, Travg = 1.0. Part (b) displays the excitation numbers, 
state synthesis fidelity and the entanglement (log- negativity) 
for the one trajectory of plot (a) that ended up in the desired 
state. 



qubits, ^ = ( 1,1,1 )^. The W-state will be generated 
with a success rate 77 given by 



Note that with the same success rate the dual W state 



|VF) = -^ (1011) + 1110) + 1101)) 



is generated (see Fig. B])). 



6 




a) 



Figure 3: (Color online) Generation of two-qubit Bell state 
|$+) = (|00) + \ll))/^/2^. (a) Probability density of the inte- 
grated (cumulative) phase shift signal X{t) — X(t')dt' 
from 6000 runs of the simulation. At time to Hadamard gates 
are applied to both qubits. Part (b) displays the probability 
density of the entanglement measure En^ the log- negativity. 
Note that neither of these plots can be obtained from the 
standard, non-stochastic master equation (see main text), i.e. 
quantum jump trajectory simulations are essential. 



Generation of GHZ-states 



Extending the 
qubits, we find 
(1,0,0,0,0,0,0,1) 



two qubit EPR scheme 
for the amplitude vector 
, and for the characteristic equation 
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Figure 4: (Color online) Generation of three qubit W-states: 
(a) Quantum trajectories for the different states that can arise 
from the given input state |^o)- At time Fto = 3, Hadamard 
gates are applied to all qubits. Windowed averaging is per- 
formed as in Fig. Part (b) displays the excitation num- 
bers, state synthesis fidelity and the log-negativity for the one 
trajectory of plot (a) that ended up in the desired W-state. 
Here trip denotes the partial trace over qubit number i, and 
the resulting pairwise entanglement happens to be the same 
for all choices of qubit pairs in this example. Note that in 
the target state all pairs of qubits are mutually entangled 
which is characteristic for the W-state and the reason for the 
robustness of its entanglement compared to the GHZ state. 



for the couplings in case of a desired GHZ state as the 



target state \GHZ) 



V2 



(|000) + |111)): 







^1 + ^2 + ^3 



which is fulfilled, for example, by the choice of coupling 
vector G = (1,-1/2, —1/2)"^. The success rate is 77 = ^. 
Again we plot phase shift signal, excitation numbers, log- 
negativity and fidelity to illustrate the correctness of our 
considerations (see Fig. ([5|). Note that due to the un- 
equal couplings, the qubit excitations and pairwise en- 
tanglement do depend on the qubit index, in contrast to 
all our previous examples, where the couplings had been 



7 



equal in magnitude. 

It is noteworthy that this 3-qubit GHZ scheme yields a 
75% chance of obtaining a W-state between qubits 2 and 
3 as a byproduct. So this might in fact be also considered 
an even more efficient scheme to generate 2 qubit W- 
states than just with two qubits in the cavity. 



EFFECTS OF DECOHERENCE 



We include decoherence into our model by considering 



the stochastic master equation Eq. (11) with the Lind- 
blad decay and dephasing rates now different from zero. 
Assuming equal rates for all the qubits, evidently entan- 
glement will be on average destroyed on a timescale set 
by T2 = (71/2 + 70)"^. When considering experimentally 
reachable parameters, which we will do further below, we 
will find that indeed the time needed to synthesize states 
is orders of magnitude shorter compared to T2 . It is thus 
clear that the simulation of the examples will look like 
above with a weak decay of coherence superimposed on 
the trajectories. 

In contrast to the decay due to decoherence, the de- 
cay due to relaxation (at a rate 71) is stochastic, in the 
sense that it leads to sudden quantum jumps. This can 
be understood by considering that the phase shift mea- 
surement stabilizes a certain subspace. Doing so, certain 
configurations of diagonal elements in the density matrix 
constitute attractors that compete with the exponential 
decay due to 71. More formally speaking, the master 
equation is unravelled with respect to the 71-process, but 
still an ensemble average description of the pure dephas- 
ing physics. The result is that the decay is stochastic 
when looking at single trajectories and the usual expo- 
nential 7i -decay is recovered when averaging over many 
trajectories. Conversely, in a single trajectory the off- 
diagonal elements decay on a timescale set by T2 , showing 
the following behavior: As long as the relaxation jump 
process has not happened, one observes a decay solely 
due to pure dephasing (see Fig. (|6|). Once the relaxation 
process has happened, coherence and thus entanglement 
are also lost completely. 
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Example - two qubit Bell-states including dissipation 

To demonstrate the influence of decoherence and re- 
laxation, we repeat the example for a two qubit Bell- 
state, 1^^) = (1 01) + 1 10)), assuming comparatively 
low values of f /71 = 10 and f = 20 to illustrate the 
effects and make all the dynamics visible. Experimental 
ratios would be at least about a factor 100 higher and 
thus the fidelity and lifetime are higher in experiment 
than they appear from the following simulations. 

We have plotted the time-evolution for the choice of 
couplings that leads to the creation of a two qubit Bell- 



Figure 5: (Color online) generation of three-qubit GHZ states: 

(a) Quantum trajectories corresponding to the different states 
that can arise from the given input state. At time to, 
Hadamard gates are applied to both qubits. Note that among 
the unwanted outcomes there are two-qubit |^^) Bell-states. 
These are actually generated with a success rate of 77 = 3/4 
which is higher than in the original two-qubit scheme. Part 

(b) displays the excitation numbers, state synthesis fidelity 
and the log- negativity for all pairs of qubits for the trace of 
part (a) that ended up in the desired GHZ-state. Note that 
once the GHZ-state is reached, all pairwise entanglement is 
lost. This is a typical feature of GHZ states, which contain 
only genuine three-particle entanglement. Part (c) shows the 
evolution for the particular trajectory that reaches the Bell- 
state between qubits 2 and 3, which can be generated very 
efficiently as a byproduct using this 3-qubit GHZ scheme. 
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Figure 6: (Color online) Effects of adding decoherence to 
the dynamics. The situation is identical to the simula- 
tion of Fig. ([2]), with the target state |^^), except for the 
added relaxation rate 71 = 0.01 • F and pure dephasing rate 
70 = 0.02 • F. We can observe that the subspace of choice is 
stabilized before the eventual decay due to relaxation. How- 
ever, even before the sudden jump due to relaxation, one ob- 
serves a slow decay of the fidelity and entanglement between 
the qubits, due to the pure dephasing rate 7^ (dashed lines). 



state |^+) = ^ (|10) + |01)). The results are shown 
in Fig. (|6|, which should be compared against Fig. ([2|. 
Likewise, we have considered the probability density for 
the time- averaged phase shift signal and the entangle- 
ment measure, for the Bell state see Fig.[7| There, 
the strict upper envelope for the entanglement is particu- 
larly noteworthy, corresponding to the decay of coherence 
within the subspace selected by the measurement. 



EFFECT OF IMPERFECTIONS DUE TO 
PARAMETER SPREAD 

In order to prepare states in this way experimentally, 
one faces the problem that it might not always be possi- 
ble to fix important parameters perfectly. If the scheme 
one has in mind in turn relies on exact matching of pa- 
rameters too much, one quickly ends up with a proposal 
that may be interesting but not very realistic. We there- 
fore examine the effects on the fidelity and entanglement 
properties of this state synthesis scheme in the presence 
of small deviations in the couplings of the qubits to the 
cavity 



9 



(16) 



where Go is the ideal coupling vector from solving the 
characteristic equations Eq. (10) for the target state. 



21 

in 



CD 



CD 

2-1 

(V 







TuT2 

N- 


10) 






00) 




01) 




preparation 







100 




dgi/g are the deviations from that ideal coupling for each 
qubit. Without loss of generality we first look at the case 



Figure 7: (Color online) Effects of adding decoherence 
and relaxation to the creation of the Bell state |$^) = 
^ (1 00) -h 1 11))- The situation is identical to the simulation 
of Fig. (jsj, except for the added relaxation rate 71 = 0.1 • f 
and pure dephasing rate 7^ = 0.05 • F. Part (a) shows the 
probability density of the time-averaged (cumulative) phase 
shift signal X{t) with an example trajectory shown in red. 
Note the build-up of finite probability at finite signal values, 
before relaxation back to zero phase shift, which represents 
the vacuum state |00) at long times. Part (b) shows the prob- 
ability density of the entanglement (log-negativity). We can 
observe that the entanglement is lost on a timescale given by 
T2 = (71/2 + 72)"^- Note in particular the sharply defined, 
exponentially decaying envelope that defines a strict upper 
bound for the entanglement at any given time. This is due to 
the pure dephasing. 
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where all couplings are equal to their ideal value, except 
one which differs by 5g/g. We focus on the stochastic 
term in the master equation Eq. (11) which is responsi- 



ble for the projection onto a set of states, one of which 
is our target state. Let us rewrite this term a bit by 
inserting the definitions of f and TV: 



Mr 



(np- 



^pN-2p(N)(t)]m 



N 



^cavity Q 



Y,^4{nip + ph,-2p{ni){t))m{^^) 



From this form of Eq. ( 17), it is evident that due to the 



linearity in the couplings we can pull out all deviating 
couplings into separate terms which have the same form. 
This reads 



N (0) 

= -V4f V^(n,p + pn, -2p(n,)(t))^(t) 

i=l ^ 

- ^ 5 

- \/4f ^^(n^p + pn, -2p(n,)(t))^(t) 

i=l ^ 

= -^/^(Np + pN-2p(N){t))m 

N 

- ^y4^(n,p + pn, -2p(n,)(t))^(t), 



which means that in addition to the ideal behaviour cap- 
tured by the first term, each individual qubit with devi- 
ating coupling will be projected on its ground or excited 
state on a timescale given by the inverse of the individual 
measurement rate 5Ti = "^^^ (we have assumed posi- 
tive 6gi for simplicity; otherwise the signs in the last line 
would change for those qubits with dgi < 0). This has 
two consequences: The first consequence concerns the 
measured phase shift: Instead of being equal for all the 
base kets that form our target state, there will be devia- 
tions in the phase shift from base ket to base ket. This 
means that we will be able to observe the breakdown of 
the target state. Therefore, second, the lifetime of the 
desired entangled state will now also be limited by the 
inverse of the maximum of the individual measurement 
rates, in addition to the effects of decoherence. 

In other words: As soon as we have gained enough 
signal to noise ratio to discriminate the different base kets 
from each other (i.e. resolve the different corresponding 
phase shifts), our target state will be destroyed. 

To illustrate this effect in a fairly drastic way, we 
choose an example of three qubits with a W-state as a 
target state and the coupling vector ^ = (1,1,1)^ + 
(72/10,0,-^2/10)^ . This yields an individual mea- 
surement rate for the second qubit of ^r2/r = 1/50. 
Therefore, we expect the target state and especially its 
entanglement properties to be destroyed on a timescale 




Figure 8: (Color online) Probability density of the time- 
averaged cumulative (integrated) phase shift signal X as an 
illustration of the effect of parameter spread in the couplings. 
Apart from the deviation in the couplings from the ideal val- 
ues, the setup is identical to the example in which we aimed 
for a three qubit W-state, as seen in Fig. |4] Hadamard gates 
are applied to all qubits at time to = 3 • f During the 
following time interval of length F"^ all trajectories are pro- 
jected onto the W-state, \W) ,|000), or |111). Meanwhile 
the competing projection on the individual number states of 
the qubits becomes more pronounced and dominates the dy- 
namics on a time scale ^F^^ = 50r~^. This is exactly the 
timescale on which we can be sure to identify all the product 
base states by their phase shift values individually. 



of 50 times the preparation time f~^. As we will ar- 
gue in the following section, present experiments allow a 
ratio T/jdecoh = O(IO^) which justifies to ignore deco- 
herence for the moment. The resulting simulation beau- 
tifully confirms the expectations, see Fig. |8] 

We conclude that in order to observe the full dynamics 
of the system one should strive for a regime where the 
condition 



F > ^F^ > 7i,70Vi e l,...,iV 



(18) 



is met. In the next section we will show that this is indeed 
possible with present-day experimental parameters. 

As a side remark we state that the situation of one 
coupling deviating from the others is in principle already 
found when synthesizing GHZ states for an odd num- 
ber of qubits (e.g. 3) as examined in the previous sec- 
tions (see Fig. (Isl a))). We had chosen a coupling vector 
G = (1,-1/2,-1/2)'. Here the lar ger magnitude of 
the coupling for the first qubit is responsible for the gen- 
eration of 2-qubit W-states. As a consequence, we can 
learn about the state of the first qubit while we can still 
not distinguish qubits 2 and 3 from each other. Follow- 
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ing our previous reasoning in this section, we find that 
the state of qubit 2 should be discerned on a timescale 
set by ((5r2/f ) = 4 (in units of f ~^), which matches 
the simulation results shown in the previous section (see 
Fig.[§i). 



POSSIBILITIES FOR EXPERIMENTAL 
REALIZATION 



Cavity QED setups in superconducting circuits [5l [47l 
HSl [49l have been used to implement ideas of quantum 
optics on the chip, and are considered a promising can- 
didate for scalable, fault tolerant quantum computing 
(e.g. [50 ). Proposals for generating and detecting non- 
classical photon states exist or have been implemented 

These experiments realize a Jaynes-Cummings cou- 
pling between qubit and resonator of up to 27r • 100 MHz, 
resonators with frequencies of about 27r • 5GHz, and a 
large range of resonator decay rates k between 10 kHz 
and 100 MHz. Given this parameter space and assuming 
a bare qubit coupling of go ~ 27r • 100 MHz^ detunings in 
the GHz range, |e|^ ~ 10 photons in the readout cavity, 
and a qubit decay rate 71 ^ 0.6 MHz^ it is easily possi- 
ble to reach values ofT/hz ^ 10^. This gives ample time 
for the state synthesis before decoherence starts playing 
a role. 

Furthermore, couplings can be adjusted with enough 
accuracy such that the state generation is also not lim- 
ited by this factor. We can examine the sensitivity 
of the ratio 6T/T to small deviations in the parame- 
ters. From ST oc Sg'^ and Sg = 6{gQ/A)^ we find 
5T/T = [26 go /go — 5A/A]'^. Assuming an uncertainty 
about the bare value of the coupling of the qubits to the 
cavity and an uncertainty about the qubit detuning of 
about 5% each we find that ST - 0.05^f . Note that this 
value is obtained without even considering the possibility 
of actively compensating for the spread in the couplings 
by suitably adjusting the detuning. This hints that un- 
der presently available optimal experimental conditions, 
the infidelity due to parameter spread becomes visible 
only long after the system has been severely decohered 
anyway. However, one can always intentionally choose 



parameters such that Eq. (18) is fulfilled and the full dy- 



namics discussed here can be experimentally observed, 
including the ultimate measurement-induced decay of the 
temporarily produced entangled state. 

The main challenging step to be taken experimentally 
before realizing this scheme in the lab is to operate in the 
single-shot qubit readout limit. This has been demon- 
trated very recently by the Saclay group using a Joseph- 
son bifurcation amplifier setup ^53j . 



CONCLUSIONS 

We have analyzed a very general experimentally di- 
rectly relevant way to generate entangled multi-qubit 
states using a dispersive phase shift measurement of the 
collective state of several qubits inside a cavity. We have 
given criteria for the possibility to synthesize a given tar- 
get state and studied the most relevant examples of Bell- 
states as well as W- and GHZ-states for two and three 
qubits. We have also discussed, and analyzed by exten- 
sive numerical simulations, the two major sources of im- 
perfections in this setup, namely decoherence and param- 
eter spread. Finally, we have compared with presently 
reachable experimental parameters and conclude that 
this scheme could soon be tested in the laboratory. 
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